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Abstract. - The simplified bandstructure introduced by Mazurenko et al to model FeSi is used to 
analyze the singlet semiconductor to ferromagnetic metal transition in the isoelectronic isostruc- 
tural alloys, FeSii-^Gea,. The complex bandstructure of the alloy is replaced by an alternating 
chain of doubly and singly degenerate atoms to represent Fe and Si/Ge respectively. The for- 
mer(latter) form narrow(broad) bands with a substantial hybridization between them. A sub- 
stantial onsite repulsion including a Hunds rule coupling is introduced on the Fe sites. The mean 
field phase diagram contains a first order phase transition from the singlet semiconductor to a 
ferromagnetic metal with increasing temperature and interaction strength similar to the alloys. 
The analysis also reproduces the rapid rise of the spin susceptibility in the semiconductor with 
a crossover to a Curie-Weiss form at higher temperatures. Good agreement is found at zero 
temperature between the mean field and accurate DMRG calculations. 
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Introduction. — The rapid crossover around room 
temperature in FeSi from semiconducting to metallic be- 
havior with strong magnetic fluctuations [l] has long 
been studied both experimentally [2j|9] and theoretically 
The small scale of the energy gap in the semi- 
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conductor has lead to proposals that this system is an un- 
usual example of a Kondo insulator in a transition metal 
^compound [2114 13 . An alternative interpretation pro- 



poses that FeSi is a nearly ferromagnetic semiconductor 
10 11 16 . A variant on this proposal interpreted the 



anomalous behavior FeSi at room temperature and zero 
field as arising from proximity to the critical point of the 
first order transition from singlet semiconductor to ferro- 
magnetic metal with increasing magnetic field. This pro- 
posal received support from a series of electronic struc- 
ture calculations using the LDA+[/ method 22 which 
found that a ferromagnetic metallic state was close by 
in energy to the paramagnetic semiconducting state 14 



Subsequently we realized that the isostructural isoelec- 
tronic compound FeGe has a magnetic metallic ground- 
with a long period rotation of the mag- 
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netic moment due to a Dzyaloshinskii-Moriya interaction. 
The isostructural isoelectronic FeSii-^Gea; alloy system 
displays an anomalous semiconductor to magnetic metal 
transition [7 18 . This transition occurs as the lattice is 



expanded. Note this behavior is opposite to the usual 
Mott transition from a paramagnetic metal to a magnetic 
insulator with increasing lattice parameter. 

The earlier electronic structure calculations are re- 
stricted to zero temperature. A phenomenological exten- 
sion to finite temperature was developed later in terms of 
interacting electrons and holes excited across the small 
band gap which gave a good description of the metal- 
insulator transition in the FeSii-zGez alloys 18 . This 
left open however the question of the underlying micro- 
scopic model responsible for this unusual behavior. It is 
not clear how to use the Kondo insulator proposal as basis 
for a microscopic model since the band structure of FeSi 
shows strong, not weak, hybridization between the Fe-3d 
band and Si valence bands. Also the magnetism of the 
metal is essentially ferromagnetic and not antiferromag- 
netic, as would be expected from a RKKY-picture. For all 
these reasons it is desirable to have a microscopic model 
closely based on the band structure which can reproduce 
the observed phase diagram. 

Simplified Multiband Model for FeSi. Al- 
though FeSi is a cubic compound it has a complex unit 
cell containing 4 formula units with low point symmetry 
at the Fe and Si sites. The ab initio LDA band structure 
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Fig. 1: (Color online) Sketch of the ID FeSi model. The dif- 
ferent line colors denote the various Fe-Fe and Fe-Si hopping 
processes. In the DMRG calculations we adopt open boundary 
conditions and terminate the chain by putting a Si orbital at 
each end. 
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calculations for FeSi correctly predict a small gap semi- 
conducting ground state 12 14.19 . This band structure 



disagrees with the Kondo insulator approach which sup- 
poses a narrow Fe 3d band weakly hybridizing with Si s, p 
bands at the Fermi energy. The semiconducting gap does 
not arise from weak hybridization between the Fe-3ci and 
Si-3p orbitals but lies with a complex of bands that are 
predominantly composed of non-bonding Fe-3d orbitals. 
The whole band structure of FeSi is quite complex so that 
it is useful to have a simpler model which contains the 
essential physics. To this end Mazurenko and collabora- 
tor 21 recently examined the Fe-Si hybridization in de- 
tail and then introduced a simplified model. Fe sites are 
represented by doubly degenerate orbitals and substan- 
tial onsite Coulomb correlations. These bands hybridize 
strongly with single orbital sites which form a broad band 
representing Si. To keep matters as simple as possible we 
analyze a one dimensional model with alternating Fe and 
Si sites. As we will discuss below, this simplification is 
not essential. It has the advantage that it allows us to 
use the powerful DMRG (density-matrix-renormalization- 
group) method [25] to check the results that we obtain 
using a mean field approximation. 

The simple multiband model is illustrated in Fig. [T] The 
Fe sites are represented by doubly red circles and Si sites 
by single dark green square. A narrow band width is cho- 
sen for the Fe-sites, a broad band width for Si sites and a 
substantial Fe-Si hybridization is introduced. A substan- 
tial Coulomb interaction with onsite Hund's rule exchange 
is included on Fe sites while interactions on Si sites are ig- 
nored. To keep matters as simple as possible we consider 
only a one dimensional model with an alternating chain of 
Fe and Si sites. Since we are interested in the competition 
between ferromagnetism and paramagnetism special fea- 
tures of one dimension such as perfect nesting are not rele- 
vant. In the mean field theory presented here we restricted 
our attention to solutions without orbital polarization, i.c 
{n m ,tj) on Fe is always independent of the orbital index m. 
Actually we found a lower energy within mean field theory 
could be achieved by relaxing this condition, e.g. allowing 
states without net spin and orbital polarizations through 
compensating changes in the spin and orbital density such 



Fig. 2: (Color online) Band structure of the simplified FeSi 
model for the case without Fe-Si hybridization (left panel) and 
with strong Fe-Si hybridization (right panel). 



that n m<a — tifn.-a^ ^m.-er = «m,cr on Fe sites. These 
states build in onsite spin-orbital correlations. We argue 
these states should not be compared to the simple mean 
field states we consider here. This conclusion is supported 
by the good agreement between the form of the mean field 
and exact DMRG phase diagrams. 

The Hamiltonian consists of kinetic energy (Hq) and 
interaction (Hi) terms. The former is parametrized by 
the hopping elements ti , . . . , £4 illustrated in Fig. [l] The 
resulting band energies £i(k) are obtained by diagonalizing 
a 3 x 3-matrix, e™" 1 , where m = 1, 2 and m = 3 denote 
Fe d- and Si p-orbitals, respectively. 
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2t$ cos k 2^4 cos k 2t\ cos | 
2i 3 cos k 2ii cos | 
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The hopping parameters from the band structure are ob- 
tained after the projection procedure 21 : ii=2.0 eV, 
i 2 =3.0 eV, t 3 =-1.2 eV, i 4 =-0.25 eV. In the following con- 
text, without special claim, all energy variables are in units 
of eV. 

The resulting energy bands are illustrated in Fig. [2j In 
the left panel the bands appear without hybridization. In 
the right panel the effect of a substantial Fe-Si hybridiza- 
tion is to split the bands into a lower bonding and an upper 
anti-bonding bands with a non-bonding band in between. 
When we choose a filling of 4 electrons per unit cell, the 
Fermi energy lies in the small energy gap between the non- 
bonding and anti-bonding bands. Thus, we obtain a small 
gap semiconducting ground state when interactions are 
ignored. Note the small value of the energy gap is not a 
consequence of a weak Fe-Si hybridization. Rather it sep- 
arates two Bloch states with predominantly Fe-character 
arising from the two orbitals on the Fe-sites. 

One consequence of the choice of a one dimensional 
model is the divergence of the density of states (DOS) 
at the band edges. We believe this is not important for 
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Fig. 3: (Color online) Zero-temperature mean field phase dia- 
gram for < U < 5eV. The dashed line (black) is the first order 
and the solid line (red) the second order transition boundary 
between the paramagnetic insulating state (PMI) and the fer- 
romagnetic metallic state (FMM). They are separated by the 
tricritical point (square blue dot). Insets: energy bands of the 
chain (a) PMI with U = 2eV, J = 0.5eV and (b) FMM with 
U = 4eV, J = leV. 



two reasons. First, as will be discussed below, the tran- 
sition between the semiconducting and the ferromagnetic 
metallic phases is a first order transition to a ferromagnet 
with a substantial moment so the exact form of the DOS 
peak is not relevant. Secondly the full FeSi band struc- 
ture calculations also show narrow peaks in the DOS at 
the band edges in the semiconducting phase. 

Mean field approximation at T=0. The inter- 
action term in the Hamiltonian is introduced only on Fe 
sites. It takes the standard multi-orbital Hubbard form 
including both intra- and inter-orbitals onsite Coulombs 
interactions 



Hi =U y^reiitniij, + n i2 t«i2j.) 
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The intra- and inter-orbital onsite Coulomb interactions 
U and U' are related through the Hund's rule exchange 
energy J by U = U' + 2 J. As we shall discuss below 
the Hund's rule exchange energy, J, plays an important 
role in stabilizing ferromagnetism. Several earlier studies 
which focussed on the anomalous T-dependence of para- 
magnetism in FeSi omitted the inter-orbital Hund's rule 
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First, we analyze the combined Hamiltonian Ho + Hi 
at T = using a mean field approximation for the case 
of 4 electrons per unit cell. The paramagnetic small gap 
semiconducting (PMI) phase remains stable initially with 
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Fig. 4: (Color online) Square of the ferromagnetic moment, 
m 2 , versus J/U at zero temperature. The tricritical point with 
J/U ~ 0.3 and U c ~ 1.7eV (blue square dot in Fig|3| marks 
the transition from a first order transition at large U (small 
J/U) to a second order transition at small U (large J/U). 



a reduced energy gap as the interactions are introduced. 
In this phase as apparent in the left panel of Fig. [3] the 
4 electrons are distributed between Fe and Si sites with 
density n Fe : risi = 3 : 1. Upon increasing the interaction 
strength there is a transition to a ferromagnetic metallic 
phase (FMM) with a moment primarily on the Fe sites, 
m Fe = X)m( n -Fe,mt - «Fc,m|)- This leads to a splitting of 
the bands between minority and majority spin directions 
as illustrated in right inset of Fig. [3| The phase boundary 
in Fig. [3] between these phases depends on the ratio J/U 
which we denote as a. Along the phase boundary between 
PMI and FMM the transition is first order at large U and 
small a but turns into a second order transition for values 
of a greater than a critical value, a c rs 0.3 as shown in 
Fig. [4] The dimensionless parameter a — J/U should 
remain almost constant passing from FeSi to FeGe. The 
value of U is to be compared to the bandwidth, therefore 
increases in FeGe relative to FeSi due to the increase in the 
lattice constant. Thus, the change between two materials 
would follow along a horizontal line in Fig. [3j The bands 
in Fig. [3] inserts, illustrate the change between a PMI and 
FMM, at PMI side in the left panel with a reduced energy 
gap for triplet excitations and the FMM side in the right 
inset panel. The transition in the magnetic momentum 
at a constant a = 0.25 and increasing U is illustrated in 
Fig. [5} Initially there is a jump followed by a rapid rise to 
a value of via — 1 which agrees with the observed moment 
of « IfiB per Fe in the magnetic phase of FeGe 7, 23, 24 . 
This unusual value is a consequence of the strong Fe-Si 
hybridization which causes only a partial polarization of 
a single band. Further increases in U > 4.1eV lead to a 
slower increase in the polarization since the Fermi energy 
in the majority spin bands enters a region with a rapidly 
rising dispersion. Note the earlier LDA+U calculations 
14 18 with the full FeSi band structure gave a similar 
moment ~ IfiB- 
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Fig. 5: (Color online) Ferromagnetic moment m versus U at 
constant J/U = 0.25 and T = 0. 



If the interaction strength is strongly increased a new 
antiferromagnetic phase is stabilized (U > 5.9eV). In this 
parameter region the strong onsite Coulomb repulsion on 
the Fe sites causes an electron to be transferred to the Si 
site resulting in the distribution np e : nsi — 2:2. The 
spins of the two electrons on the Fe site are aligned paral- 
lel according to Hund's rule. Superexchange through the 
Si sites aligns the spins on successive Fe sites antiparal- 
lel resulting in an antiferromagnetically ordered insulator. 
However, this phase which appears only at large values 
of U, is not relevant to FeSi or FeGe. Note, the relevant 
values of the onsite interaction, U, here are larger than 
the value of leV used by Mazurenko et al [21] in their 
Dynamic Mean Field Theory calculations and fits to the 
optical properties and spin susceptibility of FeSi. 

Finite Temperature Results. — We turn now to the 
temperature effects within the mean field approximation. 
The earlier model for FeSi based on LDA+U calculation 
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placed FeSi in the PMI phase with the FMM phase 
nearby at a slightly larger value for U. Extending the 
mean field results to finite T for fixed a — 0.25 and val- 
ues of U near the critical value U c = 2.6eV, leads to the 
(T,U) phase diagram shown in Fig. [6j At small values 
of T and U < U C (T = 0) there is a first order transi- 
tion from the PMI phase to the FMM phase at finite T 
driven by the larger entropy of the FMM phase. At higher 
temperatures the ferromagnetic order disappears in a sec- 
ond order transition to a paramagnetic insulating phase. 
These two transitions are illustrated in Fig. [7] which shows 
the temperature dependence of the ferromagnetic moment 
m Fe (T) at fixed a and U. As U is decreased further from 
U C (T — 0), the temperatures of the first and second order 
transitions approach each other and come together at a 
tricritical point located at T t = 0.07eV, and U t — 2.55eV. 

For small values of U (< 2.55eV), the results of calcu- 
lation of the magnetic susceptibility, x{T) are shown in 
Fig. [8j The locus of the susceptibility maxima, T max {U) 
are plotted in Fig. [6] We see that T max (U) extrapolates 
to the tricritical point. This behavior is consistent with 
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Fig. 6: (Color online) Mean field phase diagram in (T,U) 
space with J/U = 0.25 with U c ~ 2.6eV and a tricritical point 
at Ut ~ 2.55eV. The square dots (red) mark the first order 
transition between paramagnetic insulating and ferromagnetic 
metallic phase (PMI and FMM), and the empty circles (blue) 
stand for the second order transition. The filled circles (black) 
indicate the temperature of the maxima in the spin suscepti- 
bility x(T). 



explanation [14] of the maximum in x(T) observed in FeSi 
at 500K by Jacarino et al, [l] , as arising from the proxim- 
ity to a critical end point of a line of first order transition 
between PMI and FMM phase in a magnetic field. 

Density Matrix Renormalization Group simula- 
tions. — In order to test the mean field analysis we used 
a numerically exact simulation of the one-dimensional mi- 
croscopic model for FeSi defined by Eqs. ([!]) and ([2]), based 
on the density matrix renormalization group (DMRG) al- 
gorithm |25| . 

DMRG Setup — The geometry used in the computa- 
tions is displayed in Fig. [T] It consists of one Fe atom 
(containing two orbitals) and one Si atom per unit cell 
and we use a setup where the chain is terminated by a Si 
atom at both ends. 

While we are interested in the case of an average den- 
sity of 4 electrons per unit cell, one has to be careful in 
choosing the correct filling for a finite system with open 
boundary conditions, in order for the lowest two bands to 
be filled. For the chosen setup the appropriate number of 
electrons is N e = 4Y C + 2, where N c denotes the number 
of unit cells. 

Our main interest is to detect and locate the partially 
polarized conducting state reached upon increasing the 
correlations on the Fe sites. The mean-field results pre- 
dict such a phase for intermediate values of U and a larger 
than 0.2. The polarized state can be found in DMRG 
by scanning the energy in each S z sector. The degen- 
eracy of the energy with respect to S z gives the spin 
multiplet S of the groundstate. In Fig. [9] we present 
the numerical phase diagram for a system consisting of 
N c = 10 unit cells (21 atoms), where the diameter of 
the filled circle is proportional to the size of the magnetic 
moment. In addition we have measured the charge gap 
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Fig. 7: (Color online) Square of the ferromagnetization, m 2 
as function of temperature for U > U± ~ 2.55 eV with J/U = 
0.25. 
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Fig. 8: (Color online) Spin susceptibility x(T) for U < Ut with 
J/U = 0.25, the four curves are for U = 2.50, 2.51, 2.52, 2.53eV. 
The maxima of Df(T) correspond to the crossing points between 
the curves and the arrow. 
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Ajv = \[E{N e + 1) + E(N e - 1)] - E(N e ), which we take 
as a measure for insulating (An > 0) or conducting be- 
havior (Ajy = 0). The charge gap is shown for N c = 30 in 
Fig. fusing empty circles. While the charge gap at U = 
corresponds to the value of the small-gap semiconductor, 
the gap seems to close in the partially polarized region as 
expected, and then opens again and reaches large values 
for large U. The agreement of the DMRG-based boundary 
of the partially polarized phase with the mean field results 
is remarkable. 

The mean field theory predicts also an antiferromag- 
netic insulating phase in addition to the semiconducting 
and the partially polarized phases. In the one-dimensional 
model considered here, a long-range antiferromagnetically 
ordered phase is not present, but is replaced by a Haldane- 
like phase. We have checked for a specific value of the in- 
teractions (U = 10, a = 0.2) that the system has a finite 
spin gap, edge-localized spin excitations, and a two-fold 
degenerate entanglement spectrum (in the most probable 
charge sector) [26], as expected for a Haldane phase. The 
region of the phase diagram in the vicinity of the semicon- 
ductor (U — 0) has also a spin gap, however no edge spin 
excitations have been detected and the entanglement spec- 
trum does not exhibit a clear two-fold degeneracy. It will 
be an interesting question to explore in the future whether 
the two phases are adiabatically connected or separated 
by a phase transition. Both possibilities can in principle 



occur, as suggested in Rcf 27 
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In conclusion we find that the singlet semiconductor to 
ferromagnetic transition observe in FeSii-^Ge^ alloys can 
be well reproduced by a simplified one dimensional model 
of the bandstructure which retains the strong inter-atom 
hybridization of the full bandstructure. The Hunds rule 
term in the onsite interaction on Fe sites plays an impor- 
tant role in stabilizing the ferromagnetic state. The spin 
susceptibility of the paramagnetic semiconductor rises 
strongly at low temperature with a maximum at a tem- 
perature which connects on to the tricritical point of the 
Curie temperature of the ferromagnetic metal as the on- 
site interaction is increased. The reliability of the mean 
field approximation was tested against accurate DMRG at 
zero temperature with good results. 
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Fig. 9: Phase diagram obtained from density matrix renor- 
malization group simulations. Filled circles denote the size of 
the partially polarized ferromagnetic moment. Empty circles 
indicate the size of the charge gap. Two of circles have been 
labelled with the value scaling to their diameter (m = 0.5 for 
the black circles and A = 2 for empty circles). The solid lines 
show the mean field phase boundaries which are extended to 
higher values of U and include the antiferromagnetic state at 
U > 5.9eV. 
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